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Abstract 



We examine the possible states of subsystems of a system of bits 
or qubits. In the classical case (bits), this means the possible marginal 
distributions of a probability distribution on a finite number of binary 
variables; we give necessary and sufficient conditions for a set of prob- 
ability distributions on all proper subsets of the variables to be the 
marginals of a single distribution on the full set. In the quantum case 
(qubits), we consider mixed states of subsets of a set of qubits; in the 
case of three qubits, we find quantum Bell inequalities — necessary 
conditions for a set of two-qubit states to be the reduced states of a 
mixed state of three qubits. We conjecture that these conditions are 
also sufficient. 



1 Introduction 



What can we believe about some parts of a system without contradicting 
what we believe about other parts? If the system is described by a set of 
numbers, and our beliefs are the probabilities that these numbers take given 
values, then a part of the system is described by a subset of the numbers and 
our beliefs about it will be given by marginal probabilities derived from the 
probability distribution of the full set of numbers. The marginal distribu- 
tions of different parts are constrained by the fact that they all come from a 
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single set of probabilities on the full system. Bell's inequalities are an exam- 
ple of such constraints. The conclusion of the EPR argument is that a single 
quantum system like an electron has a set of numbers giving the results of 
all possible measurements, even though these cannot be measured simultane- 
ously. Wigner [23] presented Bell's theorem by considering the probabilities 
for subsets of electron observables which could be measured simultaneously 
(either directly, or by measuring the electron's partner in a singlet state), 
and showing that these subset probabilities, if they derived from a single 
probability distribution on the full set, would be constrained by inequalities 
which were not satisfied by the predictions of quantum mechanics. 

Other forms of Bell inequalities can also be understood in this way, as 
compatibility conditions on the marginal distributions of subsets. Asher 
Peres [TH] has considered this problem in complete generality, bringing out its 
formidable computational complexity. In this paper we solve the special case 
in which one is given joint probability distributions for all proper subsets of 
a set of binary variables, finding necessary and sufficient conditions for these 
distributions to be the marginals of a single distribution on the full set. 

The motivation for this study is to investigate our initial question for 
quantum systems. In this case our knowledge of the system is represented by 
a mixed state, or density matrix, and our knowledge of a part of the system 
is given by the reduced state, obtained by tracing the full density matrix over 
the rest of the system. What are the constraints on these reduced states? 
Our answer to the classical problem yields a possible answer to the quan- 
tum question, as the conditions on marginal probability distributions have 
immediate analogues for quantum states of a finite set of qubits. They can 
be translated into conditions on the density matrices of proper subsets of the 
qubits, which we prove, in the case of a system of three qubits, to be neces- 
sary for the density matrices to be the reductions of a (mixed) state of the 
full set of qubits. We conjecture that these quantum Bell- Wigner inequalities 
are also sufficient conditions. For more than three qubits the corresponding 
conditions are not even necessary; this gives rise to new separability criteria 
(the generalised reduction criteria) [23]. 

A still more general problem in classical probability, which was introduced 
by George Boole pQ, is to ask when a set of real numbers Pijk- can be 
simultaneous probabilities P(Ei & Ej & Ek & . . .) for some events Ei. This 
problem has been investigated by Pitowsky [201111], who has shown [21] that 
the problem of deciding whether the relevant conditions are satisfied is NP- 
complete. The relation to the problem considered here (and by Peres) is that 
we assume that the full sample space is a Cartesian product of finite sets and 
that our events Ei are slices of this product. 

Work on this problem appears to have concentrated on a (discrete or 
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continuous) infinity of real- valued variables, i.e. a stochastic process, in which 
case there are no conditions other than the obvious ones (see (j2.1|) below); the 
Kolmogorov-Daniell theorem ^H] asserts essentially that if these are satisfied 
for all finite subsets of the variables, then there is a stochastic process of which 
they are the finite-time marginals. The focus then is on the range of possible 
processes having these marginals. This problem for bipartite quantum states 
has been studied by Parthasarathy [T7j and Rudolph |2~2~] . 

The situation in the quantum problem for pure states is, in a sense, 
inverse to this. It is not at all easy to construct an overall pure state which 
has given marginals: there are other conditions to be satisfied in addition 
to the obvious ones [TTJ [TUJ 121 E], and there is usually only one state with 
these marginals (this is the generic situation if one is given the reduced states 
of subsets containing more than half of the total number of qubits |13| |S]). 
This can be interpreted ^Hj as meaning that irreducible n-way correlation is 
exceptional in pure n-qubit states. 

However, it is not surprising that the quantum pure-state problem should 
be different from the general classical problem, since the classical pure-state 
problem is also very different. Classically, a pure probability distribution 
consists of certainty; its marginals are also pure, the only conditions to be 
satisfied by them are the obvious compatibility conditions (j2.1j) . and the 
marginals of singleton subsets uniquely determine the overall distribution. 
The quantum analogue of the non-trivial classical problem is to ask when 
a set of subsystem states is compatible with a mixed overall state. For 
identical particles, this problem has been much studied 0, but the case of 
distinguishable particles has only recently received attention. One approach 
to it is outlined in jT3j: in this paper we suggest another line of attack. 

The paper is organised as follows. In section 2 we consider the classical 
problem and present necessary and sufficient conditions for compatibility of 
probability distributions on proper subsets of a finite set of binary variables. 
In section 3 we describe the quantum problem, prove necessary conditions 
for compatibility of reduced states of two-qubit subsystems of a system of 
three qubits, and show that the corresponding conditions are not necessarily 
satisfied for a system of more than three qubits. In an appendix we review 
other work on the quantum marginal problem. 
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2 Classical marginals 



The general classical problem is as follows. Let S = {X\, . . . ,X n } be a set 
of random variables, with Xi taking values in a finite set V{. Let A, B, ... be 
a set of subsets of {1, ... , n}, and let Sa, Sb, ■ ■ ■ C S be the corresponding 
sets of variables: Sa = {Xi : i G A}. Suppose we are given joint probability 
distributions Pa, Pb, • • • for these sets of variables. What are the conditions 
for these to be the marginal distributions of a single probability distribution 
P(xi, . . . , £„)? This means that if, for example, A = {1, . . . , r}, then 

Pa(xi, . . . ,x r ) = 22 P(xi...,x n ) 

which we write as 

P A = Xs\a(P). 

There are some obvious necessary conditions: 

Z b (Paub) = ^c(Pauc) if AHB = AnC = (2.1) 

In particular, Pa is determined by Ps if A C S. We may therefore assume 
that in our given set of subsets, none is contained in another. We will say 
that the subset distributions are equimarginal if they satisfy the conditions 
(j2.1j) . We ask what further conditions must be satisfied. 

The simplest non-trivial case — which we discuss separately, for ease 
of reading, even though it is contained in the general case which follows 
- is where S is a set of three binary variables and A,B,C are the three 
two-element subsets, so that we are considering three marginal two-variable 
distributions Pi%(x, y), P 13 (x,z) and P 2 3(y,z) where x,y,z G {0, 1}. Wigner 
|24J pointed out that these must satisfy 

Pn(x,y) < P 13 (x,z) + P 23 {y,z) (2.2) 

where z = 1 — z (but these inequalities are not satisfied by the predictions 
of quantum mechanics for the measurements of the spin components of an 
electron in three directions, where joint measurements in two different direc- 
tions are performed by measuring two electrons in a singlet state). Pitowsky 
|2"U] showed that the inequalities ()2.2|) . and the inequalities related to them 
by permuting (1,2,3), are sufficient for Pu(x, y), Pis(x, z) and ^23(2/, z) to be 
the marginals of a single three- variable distribution P(x,y,z). 

To put these inequalities in a form which has a quantum analogue, we 
regard Pi2(x,y) as a function of three variables x,y,z which is constant in 
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z, and similarly for P 13 (x, z) and P 23 (y,z). Then the functions P 12 , P13, P23 
are equimarginal if they satisfy three equations like 

P 12 (x,y,z) + P 12 (x,y,z) = Pw(x,y,z) + P 13 (x,y,z) = P^x.y.z) 

where Pi is constant in y and z. 

Now the observation of Wigner and Pitowsky can be expressed in terms 
of three- variable functions as 

Theorem 2.1. Three equimarginal two-variable functions of three binary 
variables, P\ 2 , Pi 3 and P 23 , are the two-variable marginals of a three-variable 
probability distribution if and only if 

0<A(x,y,z)<l for all x, y, z e {0, 1} (2.3) 

where 

A = 1 - P x - P 2 - P 3 + P 12 + p 13 + P 23 . 
Proof. For x G {0, 1}, define o~(x) = (— l) x , and write 

o~i(x,y,z) = a(x), a 2 (x,y, z) = a(y), cr 3 (x,y, z) = a{z). (2.4) 
Then any probability distribution P on {0, l} 3 can be written 

P = I + aai + ba 2 + ca 3 + do\o 2 + + fa 2 a 3 + go\o 2 u 3 (2.5) 

for some real constants a, . . . ,g. The marginals of P are given by 

P12 = \ + 2aa x + 2ba 2 + 2d<T 1 <7 2 , 

P 13 = \ + 20(7! + 2ca 3 + 2ea l( r 3 , (2.6) 
^23 = I + 26(12 + 2c(x 3 + 2/a 2( j 3 

and Pi = \ + 4aai, P 2 = \ + 46(T 2 , P3 = 5 + 4ca 3 . 

Hence 

A = i + 2(c/ai(x 2 + eor 1( T 3 + /a 2 or 3 ), (2.7) 

i.e. A(ar,y,ar) = P(x,y,z) +P(x,y,z). (2.8) 

It follows that the inequality ()2.3|) is a necessary condition for the existence 
of the probability distribution P(x,y,z). 

To prove that it is sufficient, note that the equimarginal condition forces 
the Pij to be of the form (|2.fij) . We have to prove that there is a value of g 
such that P defined by ()2.5|) is a positive function. Let 

Q — I + aoi + 6a 2 + ca 3 + dcricr 2 + ecria 3 + fcr 2 a 3 ; 
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then the conditions on g are 

-Q(x,y, z) < g < 1 - Q(x,y, z) if a(x)a{y)a{z) = 1, (2.9) 

and 

-1 + Q(x,y,*) <<?<Q(x,y,z) if cr(x)ff(y)<r(z) = -1. (2.10) 

But 

g = i(p 12 + p 13 + p 23 + A)-i. 

Hence the condition < A < 1, together with < P^ < 1, gives 

-l<Q(x,y,z)<l 

It follows that every lower bound is less than every upper bound in (|2.9jl for 
different values of (x, y, z); the same is true of (|2.1()jl : and every lower bound 
in (|2.10|) is less than every upper bound in (|2.9j) . 

Now suppose that a(x)a(y)a(z) = 1 and a(x')a(y')a(z') = — 1. Then in 
the equations 

a(x) = ±a(x'), a(y) = ±a(y'), a(z) = ±a(z') 

either one or three of the signs are negative. If all three are negative, then 

Q(x, y, z) + Q(x', y\ z) = \ + 2(da(x) + ea(x)a(z) + fcr(y)cr(z)) 

= A(x,y,z). 

If just one sign is negative, say the first, then 

Q(x, y, z) + Q(x', y', z) = \ + 2 (bcr(y) + ca{z) + fa(y)a(z)) 

= P2s(y,z). 

In both cases we have 

Q(x jy ,z) + Q(x\y',z')>0 (2.11) 

so that every lower bound in ()2.9jl is less than every upper bound in (|2.10j) . 
Thus there is a g satisfying all of these inequalities and giving the required 
probability distribution P(x,y,z). □ 

A classical probabilist would (probably) find it more natural to prove 
necessity from the inclusion-exclusion principle, which gives 1 — A(x,y,z) 
as the probability that X\ = x or X 2 = y or X3 = z. We have given our 
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rather clumsier proof because it connects both with the proof of sufficiency 
and with the quantum problem. 

We now move on to the general case of n binary variables x%, . . . , x n . Let 
N = {1, . . . ,n}; for subsets of N, we write A G B to mean that A is a proper 
subset of B, writing A C B when we want to allow A — B; and \A\ denotes 
the number of elements of A. 

We consider probability distributions Pa for subsets A G N, regarding Pa 
as a function of (x\, . . . , x n ) which is constant in Xi for % ^ A. If P(x\, . . . , x n ) 
is a probability distribution on all n variables, its marginal distributions Pa 
can be written in terms of operators on functions of n binary variables 
defined by 

M i f(x 1 , ...,x n ) = f(xi, ...,x n ) + f(x ± , ...,xl,.. .,x n ). 

Then 

P A = M h ...M lr P where N \ A = {i u . . . ,i r }. 
The distribution P(xi, . . . , x n ) can be expanded as 

P = Y j c A cr A (2.12) 

ACN 

where the ca are real coefficients, with C0 = 2~ n , and 

a A (xi,...,x n ) = JJo-(zi), 0-0 = 1- 

Then the corresponding expansion of the marginal Pa is 

P A = 2«-^l ^ c B a B . (2.13) 

This equation can be inverted to give cacta in terms of the marginals Pa'- 

C ^ = E 2 L lBl Pb. (2.14) 

BCA 

We can now state the generalisation of Theorem 12.11 to any number of 
variables: 

Theorem 2.2. Let Pa (A C N) be an equimarginal set of probability distri- 
butions on subsets of the variables x±, . . . ,x n . These are the marginals of a 
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single distribution P(xi, . . . ,x n ) if and only if for each subset A C N with 
an odd number of elements, 

< Y (-l) |AnS| ^B(x) < 1 (2.15) 

AUB=N 
BCN 

for all x G {0, l} n . 

Proof. To prove that the condition is necessary, suppose the distribution P 
exists and let A be a subset of N with an odd number of elements. Let 
x = (xi, . . . ,x n ) and let x' be the sequence which differs from x just in 
places belonging to A: 

~x~i if i E A 
Xi if i A. 



x L 



Then 

< P(x) + P(x') < 1. 
Expanding P as in (|2.12|) . we have 

P(x)+P(x')=2 Yl c Wx). 

|AnS| even 

Using ()2.14|1 . we can express this in terms of the probability distributions 
Pg; the result is the sum in (I2.15JL This can be verified by using (|2.13j) to 
expand (f2~T5jl : 

Y (-l) lAnBl P B = Y (-l) |AnB| 2"-l B l Y c ^d (2.16) 

N\ACBcN N\ACBcN DCB 

in which the coefficient of cd<7d is 

\A\— 1 

( _l)^B| 2 n-,B| = £ (_i r2 W-™ (^ll^oD 
BDD m=\Ar\D\ ^ ' ' ' 

N\ACBcN 

( writing m — \ A D B | ) 



r _^\ac\d\2\a\d\ |(i - ly^^i _ ^_iy A \ D i\ 
;-i)i AnD i{i-(-i)i A \ D i}, 
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so the right-hand side of (j2.16D is 

2 C£>(Td = 2 E ° D ° D 

\A\D\ odd \AC\D\ even 

since |A| is odd. Thus if the distribution P exists, the inequality ()2.15|) must 
be satisfied for each subset A with an odd number of elements. 

To show that these inequalities are sufficient for the existence of the dis- 
tribution P, we first note, as in Theorem 12. 1| that the equimarginality of the 
distributions Pa gives us coefficients cb such that 

Pa = ^2 c b°b- 

BCA 

We have to prove that the stated conditions are sufficient to ensure that there 
is a coefficient cn such that 

P = E c A (?a + c N a N 

AcN 

satisfies < P(x) < 1 for all x e {0, l} n . Writing 

Q( x ) = E c ^°"^( x )' 

AcN 

we therefore need to be able to satisfy the inequalities 

— Q(x) < cn < 1 — Q(x) whenever ctjv(x) = 1 (2-17) 

and 

-1 + Q(x') <c N < Q(x') whenever a N (x') = -1. (2.18) 
Using (j2.14j) . we can express Q(x) in terms of the distributions -Pa(x) as 



AcN BCA 



V Pb V r--nw-i B i 

nn-\B\ 1 ; 



E -P-B (_])m-\B\ fn—\B\ 

2 n-\B\ i \m-\B\ 

BCN m=\B\ V 



^ (_l)n-\B\-l 

£ 1 i-|B| P ^ ( 2 - 19 ) 



2 n-|B| 
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We will now show that the inequalities ()2.15|) imply 

-— < Q < 1 . (2.20) 

Indeed, summing these inequalities over all subsets A with an odd number 
of elements (of which there are 2 n_1 ) gives 

< Yl d s p s( x ) < 2n_1 

BcN 

where 

d B= E (- 1 ) |AnB| 

AUB=iV 
|A| odd 



(— l) r (number of s-element subsets A with \A n f?| = r 
r=o s odd and iUB = iV 



r=0 V / 



n-|S|+r odd 

1 if i? = and n is odd 
if i? = and n is even 
_ 1 )n-|B|+i 2 |B|-i otherwise 

since the sum of every other binomial coefficient in the mth row of Pascal's 
triangle is 2 m_1 if m > 1. Hence 



< 1+ ^(_l)n-|S|+l 2 |B|-lp B < 2 



which, together with (jZJDD , gives IjOljl . 

It follows from (|2.20j) that if the inequalities ()2.15|) are satisfied, then 
every lower bound is less than every upper bound in (j2.17|) . and therefore it 
is possible to satisfy all of these inequalities with a single choice of cjv; and 
the same is true of (12.1? 
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To be able to satisfy both sets of inequalities simultaneously, we need 



< Q(x) + Q(x') < 2 whenever er(x) = 1 and ct(x') = -1. 

If cr(x) = 1 and cr(x') = — 1, x and x' must differ in an odd number of places. 
Let A be the set of indices i such that Xi ^ x\\ then 0b(x) = — cr B (x') if and 
only if \A fl 5 1 is odd, so 

Q(x)+Q(x') = 2 ^ c B a B (x), 

|j4nB| even 

which, as we have already shown, is equal to the sum in (|2.15jl . Hence if 
(|2.15jl is satisfied, then Q(x) + Q(x') > 0, so no lower bound in ()2.17j) is 
greater than any upper bound in (j2.18|) : and Q(x) + Q(x') < 2, so no lower 
bound in (J2.18J) is greater than any upper bound in (|2.17|) . It follows that it 
is possible to find a suitable coefficient c^, i.e. the conditions are sufficient 
for the existence of a distribution P. □ 

The proof of this theorem suggests an alternative set of necessary and 
sufficient conditions. Define the "bit flip" operator on functions of n 
binary variables Xi G {0, 1} by 

\^if)\3'li ■ ■ ■ ? •^n) fi&l) ■ ■ ■ ) ^i—li -^ii j • • • > •^n) • (2.21) 

and for any subset A = {ii, . . . , i r }, let ka = ■ ■ • k^. Then 

Theorem 2.3. Let P A (A C N) be an equimarginal set of probability distri- 
butions on subsets of the variables x%, . . . ,x n . These are the marginals of a 
single distribution P(x\, . . . , x n ) if and only if, for all x G {0, l} n , 

-— < Q(x) < 1 (2.22) 

and, for each odd subset Ac {1, . . . , n}, 

< Q(x) + k a Q(x) < 2 (2.23) 

where 

2n-\A\ 



c _ v (-ir |Ahl P 

^ ~ 9n~\A\ A' 

AcN 



Proof. If the distribution P(xx, • • • , x n ) exists, then we can expand it in terms 
of the functions a a for subsets A C N as in ()2.12jl . and we have Q(x) = 
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P(x) - c N cr N (x), as in (l2~T31) . The inequalities (l2~T7l) and (l2~T£l) follow, 
giving 



< + Q(x') < 2 whenever cr(x) = 1 and er(x') = —1. 

This is equivalent to ()2.2Hj) . Moreover, if P exists then Theorem 12.21 holds 
and the inequalities (J2.22|) follow, as was shown in the proof of Theorem 12.21 
Conversely, the stated inequalities on Q guarantee that every left-hand 
side is less than every right-hand side in both (j2.17J) and (|2.18|) . and there- 
fore there exists a coefficient cn such that P = Q + cn&n is a probability 
distribution. As in Theorem 12.2^ the equimarginal distributions Pa can be 
expanded as 

P A (x) = ^ c b°"b( x ) 

BCA 

and then, by (l2~TT!J) . 

Q(x) = c b ctb(x) where Pa(x) = 2j c b ctb(x) 

BC7V BCA 

Hence the marginal distribution of P = Q + cn<Jn over the subset A is 
%n\a(P) = ^n\a(Q) = ^2 C ' B ° B = Pa > 

BCA 

as required. □ 



3 Quantum reduced states 

The general quantum problem concerns subsystems of a multipartite system, 
with state space H — H\® ■ ■ ■®'H ri where Hi, ... , H n are the state spaces of 
the individual parts of the system. For each subset A C N — {1, . . . , n}, we 
denote the state space of the corresponding subsystem by Ha = ®ieyi^*- 
Then the problem is: Given a set of subsets A, B, . . . and states pa, Pb, ■ ■ ■ 
(density matrices on Ha,H b , ■ ■ ■)> does there exist a state p on H whose 
reduction to Ha is pa, i-e. 

p A = tr A (p) ? (3.1) 

(Here A is the complement of A in {l,...,n}, and tr^ denotes the trace 
over Ha in the decomposition H = Ha <8> Ha-) The obvious compatibility 
conditions, corresponding to the classical conditions ()2.1)1 . are 

tr B (pAu B ) = tic(pAuc) if ADB = ADC = Q>. (3.2) 
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As in the classical case, we will call a set of states equimarginal if they satisfy 
these conditions, and we can assume that none of the subsets A, B, ... is a 
subset of any other. 

There is a further question in the quantum well as asking whether 

there is any overall state with the given subsystem states as reduced states, 
one can ask whether there is a pure state with this property. This problem 
has a simplest case for which the classical and mixed problems are trivial: if 
the given marginals are those of all one-element subsets, then one can always 
construct the classical probability distribution 

f(x 1 ,...,x n ) = . . .f n (x n ) 

with one- variable marginals fx, ■ ■ ■ f n , and one can always construct the quan- 
tum multipartite mixed state 

P = pi <g> • • • ® p n 

with one-party reduced states pi, . . . ,p n (though not for fermions: see the 
appendix). But it is not always possible to find a pure state with these 
reductions. For a set of qubits, necessary and sufficient conditions were 
found in [TTj : 

Theorem 3.1. Let pi, . . . , p„ be a set of one-qubit density matrices, and let 
Aj be the smaller eigenvalue of pi. Then there is an n-qubit pure state 
with one-qubit reduced states p±, . . . , p n if and only if Ai, . . . , A n satisfy the 
polygon inequalities 

A* (3.3) 

This result has been extended and generalised by a number of authors. 
Details are given in the appendix. 

Now let us consider the conditions for the existence of a mixed state 
with given reduced states. The simplest CclSC, clS for the classical problem, 
is a system of three qubits for which we are given two-qubit reduced states 
Pi2;Pi3;P23- The form in which we have given the classical necessary and 
sufficient conditions can be immediately translated into quantum conditions 
by replacing probability distributions by density matrices, and inequalities 
between functions (holding for all values of the variables) by inequalities 
between expectation values of operators, holding for all states — that is, 
positivity conditions on operators. We can prove that this results in necessary 
conditions for the quantum problem, and we conjecture that they are also 
sufficient. 
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We will regard the reduced density matrix of a subsystem as an operator 
on the full system by supposing that it acts as the identity on the remaining 
factors of the full tensor product state space. That is, for three qubits, we 
identify pn with pyi <g> 1, P2 with 1 <g pi ® 1, etc. Then we have 

Theorem 3.2. Quantum Bell-Wigner inequalities Suppose P12, P13, P23 
are the two-qubit reductions of a three-qubit mixed state. Then 

< (#|A|#) < 1 

for all normalised pure three-qubit states \^f), where 

A = 1 - pi - p 2 - P3 + Pl2 + Pl3 + P23- 

Proof. This can be proved in a similar way to the classical version, Theo- 
rem 1231 with the help of the antiunitary "universal NOT" operator r defined 
for one qubit by 

r(a|0) + 6|l)) = a*|l) -6*|0). 

This operator satisfies r 2 = — 1 and anticommutes with all three Pauli oper- 
ators o~i {% = 1,2,3). It is antiunitary, i.e. 

r|0) = |0), r|V) = m => ®W) = (</#)*■ (3-4) 
We extend this to three-qubit states and define 

r (5> a/3 >)|/3)| 7 ) J = (-l) a+P+7 c* aP7 \a)\m) (3-5) 

where a, (3, 7 G {0, 1} and a — 1 — a, etc. This three-qubit operator is also 
antiunitary and squares to —1, which implies the "universal NOT" property 
that it takes every pure state to an orthogonal state. It anticommutes with 
the single-qubit Pauli operators a,i <g 1 <8> 1, 1 ® <Jj ® 1 and 1 <g> 1 <g> 
Any three-qubit mixed state can be written as 

p — |l + aiOi (g 1 ® 1 + bjl ® <jj ® 1 + c k l ® 1 <g cr fe (3.6) 

+ rfijdj (g (Tj- <g 1 + eifedj (g 1 <g 0fc + /jfcl <8> (Xj <g Ojfe + ^jfc(Ti (g) Oj (g) (Xfc 

(using the summation convention for repeated indices), with real coefficients 
cti, . . . , gijk- The reduced states of p are 

Pi2 = \l + 2a { ai (g) 1 + 2fej-l <g) Oj + 2dy0-j <g> cr,-, 

P13 = ±1 + 2aj<Ti (g 1 + 2c fc l (g cr fc + 2e ifc( Ti (g cr fc , (3.7) 

P23 = |l + 2&,cr, (g 1 + 2c fc l (g cr fe + 2/j fc (7j (g a k 
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and p x = §1 + 4a i (7 i , p 2 = §1 + 4fejCr,-, p 3 = §1 + 4c fc a fc . 

Hence 

A = ~1 + 2(dijai <g> (7j <g> 1 + e^crj ® 1 <g> a 3 + ijfcl <8> (Tj ® cr*.) 
= p + r _1 p r 

since r anticommutes with single-qubit Pauli operators. Thus 

(tf |A|tf) = (tf + (*|p|*> where |1>) = r|tf) (3.8) 
> since p is a positive operator. 

Since |\1>) is orthogonal to ()3.8|) also gives 

(*|A|*) < trp = 1, 

establishing the theorem. □ 

We conjecture that the condition < A < 1 is also sufficient for the 
existence of a three-qubit state with marginals pu, pn, P23- 

In the general multipartite case, the classical compatibility conditions of 
Theorem 12.21 also have quantum analogues, namely 

0< J2 (-l) |AnB| (*|PB|^> < 1 (3-9) 

A\JB=N 
BcN 

where A C iV is a subset with an odd number of elements. However, for 
n > 3 these are not even necessary conditions for compatibility (except for 
the case A = N, n odd [23 J. The proof given above for n = 3 fails because the 
universal-NOT operator r is antilinear, not linear (which has the consequence 
that r £g> 1 does not commute with 1 <g> <Xj). We illustrate this failure with a 
counter-example for n — 4. In this case ()3.9|) becomes 

< (*|A<|*> < 1, i = 1,2,3,4 (3.10) 

where 

Al = Pl - Pl2 - Pl3 - PU + Pl23 + Pl24 + Pl34 

and A2, A3, A4 are defined similarly. But consider 

p=|$)(tf| w here = ^(|0000) + |1100)). 

We find that for this state 

Ai = |(1 ® 1 - 2P+) ® Pi ® Pi + P+ ® P ® P) 
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where P + is the two-qubit projector onto the maximally entangled state 
-j= (| 00) + 1 11)), and P and Pi are the one-qubit projectors onto |0) and |1). 

Thus Ai has a negative eigenvalue — | with eigenvector ^(|0011) + 1 1 1 1 1 ) ) . 

Since the classical inequalities are satisfied by all classical states, it is not 
surprising to find that the quantum analogues like ()3.10j) are satisfied by sep- 
arable states [2H]- Thus they constitute a set of separability criteria. These 
multipartite versions of the reduction criterion [TJl 0] have been investigated 
by Hall 

Acknowledgement We are grateful to Sam Braunstein, whose re- 
mark about the Wigner inequalities set us going in the right direction. 

A Appendix: Beyond Qubits 

Tripartite systems made up of state spaces with dimensions d% not all equal 
to 2 have been studied by Higuchi {d\ = d 2 = d%) and Bravyi {d\ = d 2 = 
2, g?3 = 4), who have found necessary conditions for a set of one-party mixed 
states to be the reductions of a pure tripartite state. Their results 
follows: 

Theorem A.l. (Higuchi 10]) Three 3x3 hermitian matrices p a (a = 1, 2, 3) 

with eigenvalues A^ < Ajj < A^ = 1 — A^ — A^ are the reduced one-qutrit 
states of a pure three-qutrit state if and only if 

a a <a> b + a c , 
Pa<a b + p c , 

la<Ctb + Pc, 
0~a <h + S c , 

e« < o~b + e c , 

(a<0~b + Cc 

and ( a < e b + r] c 

where a a = \[ a) + \ { 2 a) , p a + X^, la = \^ + A^ , 

6 a = X[ a) +2X^\ e a = 2\[ a) + \ { 2 a \ ( a = 2\ { 2 a) + \ { 3 a \ ^ = 2^ + ^ 

and {a, b, c} = {1, 2, 3} in any order. 

Theorem A. 2. (Bravyi j2]) Let p\ and p 2 be two 2x2 density matrices with 
eigenvalues X a < // — 1 — A tt (a — 1, 2), and let p 3 be a Ax A density matrix 
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with eigenvalues A 3 < P3 < v% < £3 = 1 — u± — v-i — v 3 . T/ien pi, p 2 o,nd p 3 
are the reduced states of a pure state in C 2 ® C 2 ® C 4 if and only if 

A a >A 3 + p 3 (a = 1,2), 
Ai + A 2 > 2A 3 + /i 3 + V3, 
and |Ai - A 2 | < min{z/ 3 - A 3 , £ 3 - p 3 }. 

The general version of this inequality has been found by Han, Zhang and 
Guo [9 , who, however, only proved that it is necessary: 

Theorem A. 3. (Han, Zhang and Guo 9\) Let pi, . . . , p n be the reduced one- 
particle density matrices of a pure state of a system of n particles each with 
an m- dimensional state space. Let A- a ^(i = 1, . . . ,n) be the eigenvalues of p a , 
with Xi < • • • < An"'' • Then for each pair (a, b) of distinct particles and for 
each p = 1, 2, . . . m — 1, 

p p n m—1 

i=l i=l c=l i=l 

These results can now be seen as special cases of a very general theo- 
rem due to Klyachko [TJ]. A pure state is a special case of a mixed state 
with a given spectrum (1,0, ...,0). One can consider a mixed state with 
any given spectrum and then, given a set of one-particle states, ask whether 
there is a mixed many-particle state with that spectrum which yields the 
given one-particle states. Klyachko has shown how to obtain sets of linear 
inequalities which give necessary and sufficient conditions on the one-particle 
spectra. For systems larger than four qubits, there are thousands of inequal- 
ities. Klyachko's methods belong to symplectic geometry, and are similar to 
the methods he used to solve the long-standing problem of Horn, who asked 
"What are the possible spectra of a sum of hermitian matrices with given 
spectra?" Daftuar and Hayden have used these methods to find the possi- 
ble spectra of a single reduced state pa obtained from a bipartite state Pab\ 
their paper [2j contains a very readable introduction to the relevant ideas 
from algebraic topology and symplectic geometry. There is a surprising con- 
nection to the representation theory of the symmetric group, which was also 
found by Christandl and Mitchison [4 3 ; roughly speaking, their result is that 
if the spectra of Pab,Pa and p# approximate the ratios of row lengths of 
Young diagrams A, p, u, each with boxes, then the representation of the 
symmetric group Sn labelled by A must occur in the tensor product of the 
representations labelled by p and v. 
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Finally, we note that for fermions there is a non-trivial compatibility 
problem for the one-particle reduced states of a mixed state. The solution is 
as follows: 

Theorem A. 4. (Coleman j^J) An mxm density matrix p is the reduced one- 
party state of a system of n fermions if and only if each of its eigenvalues A 
satisfies < A < \jn. 
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